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Abstract 

We compute the matter bispectrum in the presence of primordial local non-Gaussianity over a wide range 
of scales, including the very small nonlinear ones. We use the Halo Model approach, considering non- 
Gaussian corrections to the halo profiles, the halo mass function and the bias functions. We compare 
our results in the linear and mildly nonlinear scales to a large ensemble of Gaussian and non-Gaussian 
numerical simulations. We consider both squeezed and equilateral configurations, at redshift z = and 
z = 1. For z = 0, the deviations between the Halo Model and the simulations are smaller than 10% 
in the squeezed limit, both in the Gaussian and non-Gaussian cases. The Halo Model allows to make 
predictions on scales much smaller than those reached by numerical simulations. For local non-Gaussian 
initial conditions with a parameter /nl = 100, we find an enhancement of the bispectrum in the squeezed 
configuration k = k% = fc2 > k\ ~ 0.01 ft. _1 Mpc, of ~ 15% and ~ 25% on scales k ~ 1 /i _1 Mpc, at z = 
and z = 1 respectively. This is mainly due to the non-Gaussian corrections in the linear bias. Finally we 
provide a very simple expression valid for any scenario, i.e. for any choice of the halo profile, mass and 
bias functions, which allow for a fast evaluation of the bispectrum on squeezed configurations. 
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1 Introduction 



The possibility of detecting a non-Gaussian component of primordial origin in the Cosmic Microwave 
Background (CMB) temperature fluctations has been at the center of intense work during the past ten 
years, both in theoretical and in observational cosmology. Such detection would in fact provide crucial 
insight in the physics of the early Universe, possibly leading to rule-out the simplest model of canonical, 
single-field, slow- roll inflation (Bartolo et al. 2004 Komatsu et al. 2009). Current contraints from 



CMB observations are consistent with Gaussianity (Komatsu et al. 2011). Still, CMB observations with 



the Planck mission (The Planck Collaboration, 2006) will significantly improve such constraints, while 
approaching the limit for an ideal CMB experiment. At the same time, several ongoing and upcoming 
projects will map very large fractions of the matter and galaxy distributions by means of weak lensing 
and redshift measurements. This will offer the opportunity to probe and constrain, by a different mean, 



the non-Gaussianity of the initial conditions (Carbone et al. 2010; Giannantonio et al. 2011). 



In fact, for the local model of primordial non-Gaussianity (Salopek & Bond, 1990), constraints from 



the CMB and the Large-Scale Structure (LSS) from current datasets are already comparable (Slosar 



et al. 2008, Xia et al. 2011). This results from the rather unexpected effect that this specific model 



has on the bias of halos and galaxies. As shown by Dalai et al. (2008), local non-Gaussianity determines 



a significant scale-dependent correction to the linear halo bias parameter, consequence of the strong 
correlation between the primordial large-scale and small-scale matter perturbations, the latter responsible 
for the collapse of dark matter halos. On the other hand, other models of non-Gaussianity or "shapes" 
of the initial curvature bispectrum provide corrections to a bias with a lesser or no scale-dependence, 
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making the detection of this effect unfeasable by means of galaxy power spectrum measurements alone 
(Wagner & Verde, 2012 Desjacques et al. 2011a|b Scoccimarro et al. 2012 Sefusatti et al. 2012). 

However, the effects of generic non-Gaussian initial conditions on the matter and galaxy correlators, 
both in the two-point and higher-order correlation functions, are various and go beyond the halo bias 
correction mentioned above. In the first place, a non-vanishing, primordial curvature bispectrum is 
linearly evolved into a non-Gaussian component of the matter bispectrum, increasingly relevant at high- 
redshift and large scales when compared to the bispectrum induced by gravitational instability (see e.g. 



Liguori et al. 2010, and references therein). Its overall signal in future redshift surveys is large enough 



to provide constraints on non-Gaussian parameters competitive with those from the CMB ( |Scoccimarro 
et al. 2004; Sefusatti & Komatsu, 2007). Indeed, the number of Fourier modes available in large- 



volume three-dimensional redshift surveys is comparable to those corresponding to bidimensional CMB 
observations. 

In addition, primordial non-Gaussianity affects the nonlinear evolution of structure, leading to small- 
scales corrections to all matter correlators. The effect on the matter power spectrum has been explored 
computing one- loop corrections in standard cosmological Perturbation Theory (PT) in Scoccimarro (2000) 
for x 2 initial conditions and more recently in Taruya et al. (2008) for the local and equilateral models. 
It has also been studied using more recent resummation schemes in PT by Bernardeau et al. (2010) and 



Bartolo et al. (2010). Comparisons with numerical simulations can be found in Pillepich et al. (2010) 



and Smith et al. (2011). The resulting corrections are small, of the order of the percent or below for the 
allowed values of the non-Gaussian parameters. The case of the matter bispectrum has been studied in 



standard PT by Sefusatti (2009) and Sefusatti et al. (2010). Here the correction to nonlinear clustering 

100, it amounts to a few percent 
1. This is true 



adds to the linear primordial component. For a local model with /nl 
for generic triangles while reaching ~ 20% for squeezed configurations at redshift z 
also in the mildly nonlinear regime and, with small variations, at lower or higher redshift. However, PT 
techniques fail to describe the small scale regime, where matter perturbations become highly nonlinear. 

The Halo Model (HM) represents an approach to compute matter correlators in the nonlinear regime, 
beyond the reach of PT techniques. Besides, the HM easily allows for the inclusion of primordial non- 
Gaussian initial conditions. Indeed, the effects of non-Gaussianity on the matter power spectrum have 
been studied in the Halo Model framework in Fedeli & Moscardini (2010) and in Smith et al. (2011), 
where measurements from N-body simulations were compared with predictions. They find that, for a 
local model with /nl = ±100, the power spectrum is modified by about 3% at z = 1 and k ~ 1 /iMpc -1 , 
scale where the correction reaches a maximum. 

In this work we provide the first computation of the effect of primordial local non-Gaussianity on 
the matter bispectrum in the highly nonlinear regime. In particular, we focus on squeezed triangular 
configurations, where one of the wavenumbers is much smaller than the other two. Such configurations 
are relevant, for instance, in constraining non-Gaussianity from weak lensing observations. Furthermore, 
for the local model of non-Gaussianity we naturally expect the squeezed limit to be mostly affected at 
large and intermediate scales. As shown in Sefusatti et al. (2010), the effect of a local primordial non- 
Gaussianity on large scales is well described by the initial component of the matter bispectrum linearly 
extrapolated at the redshift of interest. Besides, in the mildly nonlinear regime, a local non-Gaussian 
component has a noticeable effect on the nonlinear evolution, captured in part by one-loop corrections in 
PT, an effect that is larger for squeezed configurations. We thus compare our HM predictions at linear 
and mildly nonlinear scales with measurements of the matter bispectrum in the simulations presented in 



Sefusatti et al. (2010), and find a good agreement. Moreover, the HM approach adopted here allow us to 
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reach highly non-linear scales. We also show that squeezed configurations can be captured by a significant 
simplification of the Halo Model expressions. We then provide a very simple set of formulae which allow, 
for any choice of the halo profile, mass and bias functions, for a fast computation of the bispectrum in 
the squeezed configuration. In principle, these results can be extended to describe the galaxy bispectrum 
at small scales, provided a recipe for the proper Halo Occupation Distribution. However, this approach 
is not free of certain caveats and certainly not at a satisfactory level of accuracy, so we discuss it only 
preliminarily in an appendix. 

The paper is divided as follows. We introduce the definition of primordial non-Gaussianity in Sec. [2} 
The ingredients of the Halo Model, including its non-Gaussian corrections, are discussed in detail in 
Sec.[3j In sec. 4J. we focus on squeezed triangular configurations and in Sec. |4.2| we present the mentioned 



simplified set of formulae. Equilateral configurations are discussed in Sec. 4.3 Finally, our conclusions 



are drawn in Sec. [5] We leave the description of the galaxy correlators for the appendix [A} In Appendix [B] 
we explain some technical issues regarding our calculations. 



2 Primordial non-Gaussianity 

Measurements of the CMB bispectrum have confirmed, so far, the Gaussianity of the initial perturbations 



to a high degree (Komatsu et al. , 2011 ). Thus, any departure from Gaussian initial conditions is expected 
to be small, with a possible non-Gaussian component in the curvature fluctuations roughly of the order 
of 10~ 4 or less in units of the primordial gravitational potential <F We assume that a small amount of 
non-Gaussianity can be generically described, at leading order, by a non-vanishing curvature bispectrum, 
B$(ki, &2) k^), defined by (^ki^k^i^) = (2vr) 3 i3$(A;i, k 2 , &3)£o(ki + + k^). The specific dependence 
of such bispectrum on the shape of the triangular configuration {k±, k 2 , k^} depends on the details of the 
inflationary modeQand on the post-inflationary evolution of curvature perturbations. 

In recent years, a large amount of work has been devoted to characterize the predictions of distinct 
models of inflation for the shape-dependence (or "shape" , tout court) of the curvature bispectrum (see 

r recent reviews). In this work we will focus on the local 



Chen 


, 2010 


Byrnes & Choi 


2010 



model (Salopek & Bond, 1990), which is defined in terms of a quadratic correction to the curvature 



perturbations, local in position space, given by (Komatsu & Spergel, 2001) 

*(x) = 0(x) + / NL [</> 2 (x) - <<^(x))] , 



(1) 



where </>(x) is a Gaussian random field and the nonlinear parameter /nl measures the departure from 
Gaussian initial conditions. To first order in /nl, Eq. ([!]) gives 



B$(ki,k 2 ,k 2 ) = 2 f NL P^ki) P (j> (k 2 ) + eye. 



(2) 



where P(f>{k) is the power spectrum of the Gaussian component (f>, i.e. (</>k</>k') = (2tt) 3 P<f>(k) 0£)(k + k'). 
Due the particularly simple definition of Eq. Q, this model of primordial non-Gaussianity is probably the 
most popular and best studied. For instance, it allows for a straightforward implementation of local non- 
Gaussian initial conditions in numerical simulations, performed in the last few years by several groups 



Note that a primordial Bispectrum can be generated not only from inflation, but also at preheating (Enqvist 



et al. 


20(1-] 


Chambers & Rajantie 


Figueroa et al. 


2010 
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( Grossi et al. 


2007 


Dalai et al. 


2008; 


Pillepich et al. 


2010, 


Desjacques et al. 


LoVerde & Smith 


2011 


Scoccimarro et al. 


2012 


) . Prom a theoretical point of 



well approximates the bispectrum predicted by models where a non-Gaussian component is induced by 



mechanisms acting on perturbations outside the horizon, such as the curvaton model (Enqvist & Sloth 



2002 



Vernizzi & Wands 



Lyth et al. 2003) and multiple-field inflation (Bartolo et al. 2002; Bernardeau & Uzan 



2006). 



2002 



Present CMB limits on the local /nl parameter are given by — 10 < /nl < 74 (Komatsu et al. , 2011 ). 
Since the typical size of primordial fluctuations is of the order of cp ~ 10~ 5 , it is easy to see from EqTQ 
that /nl ~ 100 corresponds to a non-Gaussian correction of the order of ~ 10 -3 . This simple observation 
justifies a treatment of non-Gaussian effects at linear order in /nLj which we will assume throughout this 
work. Since the expression of Eq. provides contributions to the curvature trispectrum of the order 
O (/nl); which we neglect, we likewise ignore any effect due to a non-vanishing initial trispectrum. The 
validity of such an approximation for |/nl| ^ 100 has been confirmed specifically for the matter and 



halo bispectrum in N-body simulations with local non-Gaussian initial conditions (Sefusatti et al. 2010 



2011) 



The most direct effect of primordial non-Gaussianity on matter correlators is given by the linear 
evolution of the primordial bispectrum. The linear matter density contrast <5k(z) is related to the (con- 
stant) curvature perturbation <3?k during matter domination via the Poisson equation, which we express 
as 5k(z) = M(k, z) introducing the function 



M(k,z)^l^- 2 T(k)k\ 
6 s L m riQ 



(3) 



Here T(k) is the matter transfer function, VL m and Hq are the matter density in critical units and the 
Hubble rate today, and D(z) is the linear growth function. 

The matter bispectrum is defined by (<5ki<5k 2 ^k 3 ) = (2ir) 3 B(ki, k 2 , /^^(ki + k 2 + k3). At large scales 
it can be approximated by its leading order (tree-level) expression in PT, given by two components as 

B{k x , k 2 , h) = Bo{k h k 2 , k 3 ) + Baih, k 2 , k 3 ) . (4) 

The initial component Bq is given by 

B (k h k 2 , h) = M(kx) M(k 2 ) M(k 3 ) B$(k!,k 2 , fo) , (5) 

where, as in what follows, we suppress the implicit time-dependence. In addition, the gravity-induced 
contribution Bq is given by 



B G (kx,k 2 ,k 3 ) = 2F 2 (k 1) k 2 )P L (& 1 )P L (fc 2 ) + 2perm. ) 



(6) 



with -Pl(^) = M 2 (k)P$(k) the linear matter power spectrum and F 2 the usual kernel representing the 
second-order solution in PT given by 



5 1 / kj k 



F 2 (k,„ k,) = - + - + ^ j (k,, • k,) + -(k, ■ k 



(7) 



We refer the reader to the review Bernardeau et al. ( 2002 ) for an introduction to cosmological PT 



5 



At smaller scales the perturbative expansion for the matter bispectrum becomes more complicated as 



at next-to-leading order several one-loop corrections have to be taken into account (Scoccimarro, 1997 



Sefusattil 2009). For these nonlinear corrections we cannot separate the effect of the initial conditions 



from the effect of nonlinear gravitational evolution as the sum of two terms as in the tree-level expression 
Eq. Q. Nevertheless the perturbative description shows that additional nonlinear corrections appear 
when initial higher-order correlators, like the bispectrum Bq, do not vanish. 

As we mentioned in the introduction, the scope of this work is to investigate how non-Gaussian 
initial conditions affect the nonlinear evolution of the matter bispectrum. While PT provides predictions 
for the nonlinear matter bispectrum with one-loop corrections improving significantly over the leading 
order terms given by Eq. ([6]), such predictions are limited to the mildly nonlinear regime (typically 



k < 0.15 /iMpc -1 for equilateral configurations at z = 1, Sefusatti et al. 2010). To go beyond it is 
necessary to resort to numerical simulations or to phenomenological approaches such as the Halo Model, 
which we will discuss in the next section. 



3 The Halo Model 



In its simplest formulation, the Halo Model (see Cooray & Sheth, 2002, for a review) assumes that all 



matter in the Universe belongs to dark matter halos, identified by their mass. Therefore, two distinct 
particles will either belong to the same halo or to two different ones. This picture allows one to compute 
the two-point correlation function (or the power spectrum) of density perturbation as the sum of two 
contributions: the 2-halo term, mainly accounting for the spatial correlations of the distribution of 
different halos in the Universe, and the 1-halo term which depends instead on the spatial distribution of 
matter inside a single halo. Clearly, while the 2-halo term is expected to describe large-scale correlations, 
the 1-halo term provides predictions in the nonlinear regime. 

More concretely, the HM expression for the matter power spectrum is given by 



P(k) = P 2h (k) + P lh (k) , 



(8) 



Smith, 2000; Ma & Fry, 2000 Scoccimarro et al., 2001) 



with the 2- and 1-halo contributions given by (Scherrer & Bertschinger 1991; Seljak, 2000 Peacock Sz 



P2h(k,z) 
Pih{k,z) 



n dmMm,z)p(k,m,z) 



P h (k,m 1 ,m 2 ) 



dmn(m, z) p 2 (k, m, z) 



(9) 
(10) 



where p is the mean matter density of the Universe, n(m) is the halo mass function with n(m)dm the 
number density of halos of mass between m and m + dm, and p(k, to, z) is the Fourier transform of the 
spatial density profile p(r, to) of a halo of mass m, 



p(k, to) = 4-7T dr r 2 p(r, to) 



sin(fcr) 



kr 



(11) 



normalized so that p(0, M) = to. The 2-halo term depends as well on the halo power spectrum, 
Ph(k, mi, m,2), describing the correlation between the centers of halos of mass ra\ and TO2- As we expect 
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halos to be tracers of the underlying matter distribution, we can assume a linear bias relation between the 
halo and the matter density constrasts, so that 5h ~ b\5. Thus, at large scales, the halo power spectrum 
can be approximated as 

Ph(k,mi,m2) = 6i(mi)6i(m 2 )PL(A:), (12) 

where b\(m) represents the linear bias function for halos of mass m. Note that for Gaussian initial 
conditions, b\ only depends on the mass m, as implicitely assumed in Eq. (12). However, if primordial 
non-Gaussianity is considered, b\ will also depend in general on the scale k. Either way, the 2-halo term 
can be rewritten as 



P2h(k,z) 



1 

V 2 



.1=1 



J~J / dm,in(mi, z) p(k,m, z) b\(m, z) 



PUk) 



(13) 



with an additional dependence on k in the b\ function for non-Gaussian initial conditions. 

This description can be easily extended to the matter bispectrum. In the case of a three-point 
function, we should account for the possibility that the three points belong to just one, two or three dark 
matter halos. This means that there are now three distinct contributions to the Halo Model expression 
for the matter bispectrum, that is 



B(k x , k 2 , k 3 ) = B 3h (ki, k 2 , k 3 ) + B 2h (k 1 , k 2 , k 3 ) + B lh (ki, k 2 , k 3 ) , 



where 



B3h{ki,k2,k 3 ,z) 
B 2h (ki,k 2 ,k 3 ,z) 

Bih(h,k 2 ,k 3 ,z) 



i=l 



B h (ki,mi; k 2 ,m 2 ; k 3 ,m 3 ; z) 



dmn(m, z) p(m, z, ki) dm' n(m , z) p(m ', z, k 2 ) p(m ', z, k 3 



x P h (ki,m,m',z) + eye. , 

[dmn(m, z) p(ki,m, z) p(k 2 , m, z) p(k 3 , m, z) 
pi j 



(14) 

(15) 

(16) 
(17) 



In this case, while the 2-halo term depends on the halo power spectrum as in the previous case, the 3-halo 
term involves the halo bispectrum, Bh(k\,m\ \ k 2 ,m 2 ; k 3 , m 3 ; z). Assuming a local bias relation between 
halos and matter, 8h(rn) = f(S), expanded perturbatively as 8h{rn) = bi{m)5 + [b 2 (m)/2]5 2 + 0(5 3 ), it is 
possible to derive the tree- level expression for the halo bispectrum, valid only in the large-scale limit, in 



terms of the matter power spectrum P(k) and bispectrum B(k\, k 2 , k 3 ). This reads (Fry Sz Gaztahaga 
B h ( ki , mi ; k 2 , m 2 ; k 3 , m 3 ; z) 



1993) 



61 (mi) 61 (m 2 ) h(m 3 ) B(ki, k 2 , k 3 ) 
+ [61 (mi) 61 (m 2 ) b 2 (m 3 ) P(h) P{k 2 ) + eye] 



(18) 



where b 2 (m) is the quadratic bias function. For Gaussian initial conditions 61 and b 2 are scale independent. 
Moreover, since this equation is valid on large scales, we can replace the matter power spectrum P by its 
linear prediction Pl and the matter bispectrum B by its gravitational contribution Bq, Eq. ([6]), assuming 
a vanishing initial bispectrum Bq = 0. If the initial conditions are not Gaussian the matter bispectrum 
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B also contains a non- vanishing initial contribution Bq, Eq. §5§, and the expressions of the halo power 
spectrum Eq. ( 12 ) and of the halo bispectrum Eq. ( 18 ) are modified, particularly for local non-Gaussianity, 
to account for scale-dependent corrections to both the linear and quadratic bias functions. 

Halo Model predictions for the real and redshift space matter bispectrum with Gaussian initial con- 
ditions have been recently studied and compared to numerical simulations by Smith et al. (2008) and 



Valageas Sz Nishimichi (2011). Smith et al. (2008) presents a more detailed treatment of the dependence 



of halo model contributions on the large-scale PT expressions for the matter correlators which involve an 
explicit smoothing of the matter density by some filtering function. Here, for simplicity, we will ignore 
this issue as it does not affect the target of our calculations, that is the small-scale corrections induced 
by primordial non-Gaussianity. The model proposed by Valageas 8z Nishimichi (2011) differs from our 
as it is based on a Lagrangian description of halo correlations. Such description presents the advantage 
of predicting a 1-halo contribution that vanishes in the large-scale limit, as opposed to the constant 
asymptotic behavior resulting from Eq. (17), which has no physical justification. In fact, this large-scale 
1-halo contribution is also present in the matter power spectrum and it is a well recognized problem 



of the halo model, for which no clear solution has been proposed so far (Cooray & Sheth, 2002 Smith 



et al. 2003; Crocce Sz Scoccimarro 2008). Since here we are mainly interested in small scales effects 



we do not address it in our expressions. Finally, another issue related to the halo model description is 
given by halo-exclusion, that is the explicit accounting of the finite size of dark matter halos in halo 
correlation. For instance, a halo-exclusion correction has been implemented by Smith et al. (2011) for 
the matter power spectrum, resulting in a few percent correction on intermediate scales (k ~ 1 /iMpc -1 ). 
In that reference it is argued that the same correction could also approximately cancel the large-scale 
contribution from the 1-halo term. We have attempted to implement a similar correction for the halo 
bispectrum, but the resulting expression does not allow for an efficient numerical evaluation and we do 
not include it in our predictions. 

It should be noted that many of the theoretical improvements proposed in the literature to the basic 
halo model description of Eq. ( 15 )-( 17) ultimately fail to reduce the discrepancy of the predictions with 
simulations results below a 10% level in the mildly nonlinear regime. We therefore limit ourselves to the 
simplest expressions and we leave any refinement for future work. The final results of this paper will in 
fact consist in small scales effects of non-Gaussian initial conditions, whose description already present a 
significant degree of complexity. Aside for the large-scale approximations for the matter correlators, we 
can clearly identify three essential ingredients in the Halo Model prescription for nonlinear correlators: 
the halo mass function, the halo bias functions and the halo profile. In the rest of this section we will 
describe each of them individually paying specific attention to the corrections induced by non-Gaussian 
initial conditions. 



3.1 Halo Profiles 

In the Halo Model formulation assumed here, each halo is labeled by a single parameter: its mass. 
The spatial distribution of matter in a halo of mass m is specified by the halo density profile p(r,m), 
interpreted as an average over all halos of the same mass. 



We will consider the Navarro, Frenk & White (NFW) form for the halo density profile (Navarro et al. 



1997), 



P{r) (r/r 8 )(l + r/r a )*' (19) 
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which assumes a universal profile as a function of r. The parameters r s and p s can be expressed in terms 
of the virial mass of the halo m and the concentration parameter c. In particular, the virial mass is 
given by m = (47r/3) -R 3 p, with R v the virial radius, defined as the radius of a sphere within which 
the mean density of the halo is A„ times that of the Universe, we take A v = 200. The concentration 
parameter c is defined as c = R v /r s and is typically a function of m. Here for c we assume the form of 



Bullock et al. (2001). Thus 



3m 



4tt (■''' A r (> 



1/3 



Ps 



-A v pc 3 pn(l + c)-c/(l + c)]" 



(20) 



where the second relation has been obtained by integrating Eq. (19) over the volume and using the 
definition of m above. 

So far, the effects of primordial non-Gaussianity on the profiles of dark matter halos have not been 
extensively investigated. To the best of our knowledge, the only measurements of halo profiles in numerical 



simulations with non-Gaussian initial conditions are those presented by Smith et al. (2011). In the case 
of a local /nl = ±100, they have found that the ratio between the non-Gaussian to Gaussian halo profile 
in real space can be fitted by 



PNc(m,r,z,f N L ) 
p G (m,r,z) 



oc 



1 + y(m) log 10 



r x {m) 



(21) 



The mass-dependence of the functions y(m) and rx(m) is obtained by fitting the ratio above over the 
limited mass-range ~ 10 13 — 10 15 M Q , and further imposing y(m m i n ) = rx{m m \ n ) = 0, where m m ; n is 
the minimal mass, typically of order m m i n ~ (10 4 — 10 6 ) M , adopted in the numerical integrations 
in the Halo Model terms^j Note, however, that this constraint has no real physical motivation. 

We have checked that both functions y(m) and rx(m) approximately scale as oc R = [3m/(47r/o)] 1//3 , 
as one would expect. In particular, the data points in Fig. 8 of Smith et al. (2011) can be fitted by the 
linear functions y = sign(/NL)ai/-R and rx 
a x = -0.00188 for / NL = 



a y = 0.110 



a x R with a y = 0.136, a x = 0.00165 for / NL = +100 and 
100. From this observation and assuming that the effect of local 



non-Gaussianity is proportional to /nl> we decide to parameterize the ratio (21) by the following ansatz 



PNG{m,r,z,f NL ) 
p G (m,r,z) 



oc 



/nl 
100 



a y R(m) log 10 



axR(m) 



(22) 



where the proportionality constant is choosen to ensure that the Fourier transform of the corrected profile 
verifies p-^G(m,k — > 0, z,/nl) - > m - Only more accurate simulations will be able to confirm or exclude 
this ansatz with more precision^ The left panel of Figure [l] shows the data points corresponding to 
the quantities a y and ax derived from Fig. 8 of Smith et al.| ( |2011 ) , plotted as a function of the halo 
mass and compared to our constant, best-fit values. For such quantities, no dependence on mass is 
required by the data. On the right panel of Figure [T] we show the ratios of non-Gaussian to Gaussian 



2 This artificial low-mass cut-off must be introduced to avoid the computationally challenging integration of 
pushing the lower limit down to zero in the halo term expressions. See Appendix [B| for details. 

3 While the NFW form for the Gaussian profile allows for an analytic Fourier transform, the non-Gaussian 



correction of Eq. (21) or Eq. (22) does not have a closed expression in Fourier space. We are then forced to obtain 



p(k,m, z, /nl) numerically. 
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Figure 1: Left panel: measurements of the parameter a y = y{m)/R {top) and of the pivot point ax = 
rx(Tn)/R (bottom) describing the non-Gaussian correction to the halo profile, Eq. (21), derived from the 
plots of y(m) and rx(jn) shown in Fig. 8 of Smith et al. (2011 ), as a function of the mass m. Red and blue 
data points correspond respectively to /nl = +100 and —100 with the red, continuous lines and blue, 
dashed lines being, respectively, the constant interpolating values. Right panel: Gaussian halo profile 
normalized to unity at large scales, p(k,m)/m (top), for m = 10 12 , 10 13 , 10 14 , 10 15 and 10 16 /i -1 Mq 
(right to left) and the ratio between the non-Gaussian and Gaussian halo profiles (bottom) for the same 
masses (larger mass implies a larger correction) with /nl = ±100. All curves are obtained for z = 0. 



halo profiles in Fourier space. They correspond to redshift z = and masses ranging from m = 10 to 
10 16 h~ l M , and were computed for /nl = ±100. As expected, the transition from p(m,k)/m = 1 to 
p(m, k)/m ex l/(kr s ) 3 , around k ~ l/^s, is not modified in the presence of non-Gaussianity. Comparing 
the profiles with /nl = —100,0,100, for the same mass, we see that their relative difference is only 
noticeable in the large-momentum tails. The shift in amplitude is more appreciable for large masses. 

We can anticipate that the effect on the matter bispectrum due to the non-Gaussian corrections to the 
halo profiles is quite subdominant, particularly when compared to the effect of non-Gaussian corrections 
to the other ingredients of the HM, which we discuss in the rest of this section. 



3.2 Halo Mass Function 



The halo mass function n(m) characterizes the number density of halos per unit mass. Theoretical 
predictions for this quantity go back to the early ansatz of Press & Schechterl (1974) and have been 



largely improved over the years, either by comparing with N-body simulations (e.g. Jenkins et al. 1998 



Warren et al. 


2006, 


Crocce et al. 


2010 


Courtin et al. 


2011 


) and by more rigorous theoretical modeling 



in the framework of the excursion set theory (Bond et al. 1991 Sheth & Tormen 1999; Maggiore & 



Riotto, 2Q10a|[b Corasaniti &; Achitouv, 2011b|a Paranjape et al. 2012). 

The fraction of the total mass of the Universe contained in all the halos with mass in the range 
[m, m + dm] can be written as 

— n(rn)mdrn = f(v)dv. (23) 
The function f{v) presents an approximately universal form and depends on the variable 



5c 
a(m) 



(24) 



with 5 C representing the critical density for spherical collapses (we assume the Einstein-de Sitter value 
5 C = 1.68) while a(M) is the r.m.s. of matter fluctuations in spheres of radius R = (3m / 'Air p) l/ ' 3 (associated 
to the Fourier transform Wn(k) of the top-hat function in real space), 



a 2 (m) = Airjdk k 2 P L (k) W R (k) . 



(25) 



Since we will not be able to compare our results with numerical measurements on very small scales, we 
will omit implementing the most recent theoretical prescriptions for the mass function. Instead, we will 
adopt the Sheth & Tormen (ST, Sheth & Tormen, 1999) expression, i.e. 



f(y) = A 



a V 
~2tT 



1 



1 



(a v 2 



-av 2 /2 



(26) 



where a = 0.707 and p = 0.3 while A = 0.322 ensures a proper normalization. 

The expression above will be modified in the case of non-Gaussian initial conditions. In particular, 
since primordial non-Gaussianity affects the high-mass tail of the halo mass function, such effect has 
been early recognized as a possible test of the initial conditions (see Sefusatti et al. 2007, for a brief 
introduction to the early literature). Moreover, thanks to several new results from numerical simulations 



assuming local non-Gaussian initial conditions (Dalai et al. 


2008 


Grossi et al. 


2009, 


Pillepich et al. 


2010, Desjacques et al. 2009; Wagner et al. 2010), this topic has been recently subject of an intense 


theoretical activity 


(Lo Verde et al, 2008| 


Afshordi & Tolley 


2008; 


Lam &; Sheth 


20091 |Valageas[ 12010) 


Maggiore k, Riotto 


2010c||De Simone et al.H2011| D'Amico et al. 2011[ Paranjape et al. 


|2011||Achitouv 


& Corasaniti 2012 


Musso & Paranjape 


2012). 



We parameterize the effect of primordial non-Gaussianity in terms of the ratio between the uq and 
tinG) respectively the mass function for Gaussian and non-Gaussian initial conditions, 



Rng(iti, z,/nl) 



n NG (m,z,f] 



NLJ 



n G (m,z) 



(27) 



Moreover, we assume for such function the expression proposed by Lo Verde et al. (2008), based on an 



Edgeworth expansion about the Press-Schechter mass function. At linear order in /nl, this is given by 



#ng(^,/nl) 



1 + -u(u 2 
6 V 



3) ssH 



1 



dInT + 0{/NL) ' 



(28) 
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where S3(V) represents the reduced skewness of the initial, smoothed, matter density field £r(x) 
as 

1 'lil 



(4) 3/s 



(m) 



defined 



(29) 



with the halo mass m and radius R related to the variable v by Eq. (24). The third order moment (5^) 
is computed 4 ! as an integral over the initial bispectrum, Eq. (15]). In addition, to improve the agreement 



between the measurements of the mass function correction in the numerical simulations of Desjacques 



et al. (2009) with the prediction of Eq. (27), we introduce a scaling parameter q defined by -Rng(^) 



Rng(qv) with q = 0.91 (Sefusatti et al. 2011). 



3.3 Halo Bias 



The Halo Model expressions for the matter power spectrum and bispectrum involve, at large scales, a 
description of halos correlators, related in turn to the linear matter power spectrum and the tree- level 
matter bispectrum by an expansion of the halo bias relation. 



In the context of Gaussian initial conditions, the peak-background split approach ( |Cole fc Kaiser 



1989, Mo & White, 1996) provides a relation between the local halo number density and the large-scale 



matter perturbations, which can be turned into a dependence of the halo density contrast, Sh, on the 
local value of the matter density contrast, 5, smoothed on a scale larger than the typical size of the halos 



considered. In this case, the relation 8h(5) can be expanded in a Taylor series as (Fry & Gaztahaga 



1993) 



5 h {5) = h(m) S + - b 2 (m) S 2 + . . . , 



(31) 



which defines the linear 6i(m) and nonlinear bi>2{m) halo bias functions. The bias functions can be 
derived from the unconditional halo mass function and, in the case of the ST form, one obtains for the 



first two the expressions (Sheth & Tormen, 1999 Scoccimarro et al. 2001) 



b x (y) 



1 + 



a v 



1 



Sr. 



+ 



2p 



21 



[61(1/) -1] + 



5 c (l + (a^)P) 
3 



av 2 av 2 



5c 



1 + 2p av 



Sr. 



2p/S c 



l + (av 2 )P ' 



(32) 
(33) 



where v = v{m) is defined by Eq. (24) while a = 0.707 and p = 0.3 are the ST values of two mass function 
parameters. The requirement for the total matter density to be given by 



p [1 + <5(x)] = / dmmn{m) [1 + 5h(rn)\ 
dmmn(m) 



4 In the evaluation of the mass function correction we make use of the following fit, 

s 3 (jy) = / NL exp [-7.98 - 0.177 In v + 0.0186 In 2 v - 0.00260 In 3 v] 



(34) 



(30) 



noticing that the dependence on scales - i.e. on the variable v - is relatively weak over the relevant range. 
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imposes the condition 



Jdmmn(m) = p, 
along with the constraints on the bias functions, 



dm m n(m) b\(m) = Jdv f(y)b\{y) = 1 , 
-Jdmmn(m)bi(m) = Jdv f(v)bi(v) = 0, V i > 1 



A 
1 

A 

where <5y is the Kronecker delta. Such relations assure that, on large scales {k — > 0, p 



m 



(35) 

(36) 
(37) 
the 2-halo 



term of the power spectrum, Eq. ([9]), reduces to the linear power spectrum and the 3- halo term of the 
bispectrum, Eq. (15), reduces to the large-scale matter bispectrum. 



The expansion (31 ) allows to derive perturbative bias expansions for halo correlation functions. Equa- 



tions (12) and (18) represent the leading expressions for the halo power spectrum and bispectrum, re- 



spectively. These predictions change significantly for non-Gaussian initial conditions. Indeed, since the 
bias functions are derivatives of the mass function, any correction to the latter induces a corresponding, 
scale-independent, correction to the bias. 

As first pointed-out by Dalai et al. (2008), an initial local bispectrum is responsible as well for a 



scale- dependent correction to the linear bias. This rather unexpected effect received in the last few years 
a lot of attention since it allows to place tight constraints on the local /nl parameter from galaxy power 
spectrum measurements in current data-sets (see, e.g. Slosar et al. 2008). Several descriptions of the 



halo bias correction due to local non-Gaussianity have been proposed in the literature (see the review in 



Desjacques Sz Seljak, 2010, and references therein). Here we will follow the approach of Giannantonio & 



m 



Porciani (2010), because it allows to easily derive an expression for the halo bispectrum (Baldauf et al. 



2011 ), whose validity at large scales has been recently confirmed by comparison to numerical simulations 



Sefusatti et al. (2011). In this description, for local non-Gaussian initial conditions, the expansion 
(|31[) is replaced by a bivariate bias relation where the halo density contrast 5h depends not only on the 
matter overdensity 5, but also on the Gaussian component of the curvature perturbations 0, defined by 
Eq. ([!]). At linear order in /nl and up to second-order terms in 5 and (j), we have 



5 h (5, (p) = b 10 5 + b i 4> + - b 20 5 



+ 



&h<$0 + 0(/Sl) 



(38) 



where the bias parameters bij, with the indices i and j corresponding to the powers of 5 and 4>, respectively, 
depend on the halo mass. From the bias relation above we can derive a rather lengthy expression for 
the halo bispectrum (Baldauf et al. 2011 Sefusatti et al. 2011). This expression can be considerably 



simplified by replacing the Gaussian component of the curvature perturbation (j) by its full non-Gaussian 
value This is well justified for the linear term boi<j), while it is an approximation for the term bi\5 <fi (see 

Note however that, as &oi and b\\ are proportional 
and thus consistent with our treatment. 



Scoccimarro et al. 



2012, for a detailed discussion). 

nl) 



to /nL; this approximation is correct up to 0(f£ 

Thus, we consider the following bias relation in Fourier space 



5 h (m,k) 



h(m, k) 5 k + i J d 3 qi d 3 q 2 S D (k - q 12 ) 6 2 (m, q%, q 2 ) 5 qi <5 qi + C(f\ 



nl) 



(39) 
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where the non-Gaussian bias functions are now scale-dependent. In particular, we have 



61 (m, k) = b lyG {m) + AfeStcM + A&i° NG (m, k) + 0(f } 



(6) 



nl; 1 



(40) 



where b\ yG {m) corresponds to the Gaussian value given by Eq. (32), while the second and third terms on 



the r.h.s. represent, respectively, a scale-independent and a scale-dependent correction (Desjacques et al. 

v d 



2009 Dalai et al. 2008), given by 



Ab 



(a) 
1,NG 



(m) 



Ab 



(b) 

1,NG 



(m, k) 



5 C dv 
2/nl5 c (6i,g 



lni? NG (^, /nl) , 
1) 



M(k,z) 



Similarly, for the quadratic bias function b2(m,ki,k 2 ) we have 



62 (m, fei, fe 2 ) = b 2iG (m) + a4°ngM + A 4 b NG( m ' + C(/nl) > 



(41) 
(42) 

(43) 



where b 2yG (m) is the constant quadratic bias for Gaussian initial conditions, Eq. (33), while the non- 



Gaussian correction is given again by a scale-independent and a scale-dependent contribution ( Sefusatti 



et al. 2011 Giannantonio & Porciani, 2010), 



Ah 



(a) 
2,NG 



(m) 



1 



d 2 R^ G (iy,f NL ) 



51 #ng(^,/nl) 



A6f NG (m, fci, fc 2 ) 



2 /nl <5 c 



^2,G + 



cV 2 
13 1 
21 ~~ J r 



+ 2 (61 G - 1) Ab 



G 



1) 



1 



(a) 

1,NG ' 



+ 



M(ki,z) M(k 2 ,z) 



(44) 
(45) 



Note that the conditions (35)-(37) are satisfied also when taking into account the corrections to the halo 

tave 

(46) 



bias functions and the correction to the mass function, Eq. (27). In particular, for /nl / we still have 

' dmmn^ G {m, z, /nl) = P, 



and 



As noted by Smith et al. (2011), the last relation implies 

1 
P 



dm m nNG 



involving only the scale-independent corrections, and also 



dm m n NG A6^ NG 



(47) 



(48) 



(49) 



involving, instead, the scale- dependent part. At leading order in /nl, Eq. (49) is automatically satisfied 
by the corrections in Eqs. (42) and (45), which are proportional either to [bi jG (m) — 1] or to b 2yG (fn]R 
We have now all the ingredients to compute the matter bispectrum for local non-Gaussianity. 



Smith et al. 



(2011) point-out that, in general, the condition (49) is satisfied, for i = 1, if A6^ G (m, k) 



(&i,g + Ab\ ^ G — 1). Taking into account this and the fact that at large k the correction is negligible, they set 

A&i ng — m their calculation. While this approach is well justified, we keep the scale-dependent corrections to 
both linear and quadratic bias in our evaluation. 
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4 The Matter Bispectrum 



Let us now use Eqs. (14)- p7[) to compute the matter bispectrum. Including the corrections due to local 
non-Gaussianity, from Eq. (|39|) we can derive an expression for the non-Gaussian halo bispectrum in 
the tree-level approximation and plug it into Eq. (15). This is given by Eq. (18) with the bias functions 
&i(m) and b<i (to) replaced by their non-Gaussian scale-dependent versions given in Eqs. ( |40| ) and (43), 
respectively. Furthermore, at the tree level the halo power spectrum to plug into Eq. (|16l) 



is given by 



Eq. (12), where for the linear bias we use Eq. (40). 



In the following, all mass integrals are computed numerically. Due to the divergent behavior of the 
integrands for to — > 0, they are computationally expensive. For this reason, we generalize a numerical 



prescription used to compute the matter power spectrum with the HM ( Refregier & Teyssier 2002 ; Fedeli 



Sz Moscardini, 2010). We discuss this prescription in Appendix \B\ 



Since we compare our predictions with measurements of the matter bispectrum in the numerical 



simulations presented in Sefusatti et al. (2010), we assume throughout the paper the same cosmology 
as these simulations: a flat ACDM model with matter and baryon content given by O m = 0.279 and 
Qb = 0.0462, a spectral index n s = 0.96 and a r.m.s. of fluctuations in spheres of 8 Mpc of as = 0.81. 
We compute our predictions up to k = 10 h Mpc -1 , at redshifts z = and 1. However, the bispectrum 
measurements that we consider are limited to larger scales, i.e. k < 0.3/iMpc -1 . 



4.1 The squeezed limit 

Fig. [2] shows the matter bispectrum for isosceles triangles, B(k±, k, k), with fixed k\ = 0.014 /iMpc -1 , as 
a function of k. Thus, this subset corresponds to increasingly squeezed configurations, starting with an 
equilateral one. Left panels show the results at redshift z = 0, while right panels at z = 1. Dotted, dashed 
and dot-dashed curves correspond, respectively, to the 1-halo, 2-halo and 3-halo contributions, while 
continuous lines correspond to their sum. In the first row we present the matter bispectrum for Gaussian 
initial conditions, together with the ratio between the measurements in N-body simulations and the model. 
The shaded area corresponds to a 10% agreement. Our model is tested against simulations only in the 
mildly nonlinear regime, where we observe agreement within a 10% error. The second and last rows show 
the correction due to non-Gaussian initial conditions. In particular, the second row shows the difference 
between the predictions with and without non-Gaussianity, A£?ng = -Bng(/nl = 100) — -Bg(/nl = 0), 
while the third row shows their ratio, -Bng(/nl = 100)/I?g(/nl = 0). As shown in the middle row, the 
ratio between the N-body simulations and our predictions is inside the shaded area corresponding to a 
10% agreement for z = and slightly outside for z = 1. This is a very good agreement, if we consider 
that the non-Gaussian correction for /nl = 100 is per se of the order of ~ 10% (~ 20%) for z = (z = 1). 

At large k the bispectrum is dominated by B2h{k\, k, k). This is not surprising because this term is 
dominated by the position-space configuration in which two points are in the same halo and the third 
one is far away. We will take advantage of this fact in the following section. 

Notice also that the relatively large corrections at intermediate scales due to both non-Gaussianity 



and nonlinear evolution, studied by Sefusatti et al. (2010), extend according to the Halo Model predictions 
well into the nonlinear regime, k > 0.3 h Mpc - . Similarly to what happens for the effect on the matter 



power spectrum (Fedeli & Moscardini 2010 Smith et al. 2011), the non-Gaussian to Gaussian ratio 

1 h Mpc - . In particular, such maximum represents a distorsion 
0.8 h Mpc -1 , for z = 0, and of the order of ~ 25% at 



reaches a local maximum around k 
of the Bispectrum of the order ~ 15% at k 
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Figure 2: Squeezed configurations of the matter bispectrum, B{k±, k, k), with k\ «C &2 = k$ = at 
redshift z = {left panels) and z = 1 {right panels). See text for explanation. 
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k 1.4 /i Mpc , for 2 = 1. The enhancement in the matter bispectrum (~ 15% — 25%) due to the 
presence of local non-Gaussianity is therefore almost an order of magnitude greater than the analogous 
enhacement in the power spectrum (~ 2.5% — 3.5%). 

In order to understand the role played by each of the non-Gaussian corrections affecting the Halo 
Model ingredients, in Fig. [3] we show such individual components both for the matter power spectrum (left 
panel) and for the squeezed configurations of the matter bispectrum (right panel). For the matter power 
spectrum the corrections to the halo mass function dominate the intermediate scales k ~ 1/iMpc -1 , 
enhancing the power on these scales and inducing a bump in the ratio between the non-Gaussian and 



Gaussian power spectrum (Fedeli & Moscardini, 2010 Smith et al. 2011), mentioned above. The correc- 



tions in the halo profile are only important on smaller scales, while those to the linear bias are always 
subdominalj^} For the matter bispectrum, the initial non-Gaussian conditions dominate its amplitude at 
large scales, as one can see on the right panels of Fig. [3} However in this case, as opposed to the power 
spectrum, the corrections due to the linear bias function b\ dominate the intermediate and small scales, 
i.e. k > 0.1 /i Mpc" 1 . At k ~ 1/iMpc" 1 , also the corrections to halo mass function become important. 
Remarkably, over the range of scales considered here, the correction on the density profile is subdominant 
with respect to the corrections to the mass function or halo bias. 

In conclusion, the enhancement in the power spectrum due to the non-Gaussian corrections, is also 
present in the bispectrum. However, in the bispectrum it is dominated by the non-Gaussian corrections 
to the linear bias, and not by the corrections to the mass function as in the power spectrum. This trans- 
lates into a bump in the ratio between the non-Gaussian and Gaussian bispectra, reaching a maximum 
amplitude in the squeezed configuration around k ~ 1 /iMpc -1 of ~ 1.15 for z = and ~ 1.25 for z = 1. 

4.2 Approximating the squeezed limit 

The squeezed limit allows for a significant simplification of the Halo Model expressions. As already 
mentioned, in this limit the largest contribution to the bispectrum comes from the position-space config- 
uration where two points are close and belong to the same halo while the third one is at larger distance 
from the first two, and hence is likely to belong to another halo. In this case we expect the Halo Model 
prediction to be dominated by the 2-halo contribution, with the 1-halo and 3-halo terms being subdomi- 
nant. In Fig.jgJ where fa = 0.014 h Mpc -1 , it can be clearly appreciated that -B2/1 becomes dominant over 
the other two terms for k > 0.3 h Mpc" 1 . At k > 1/iMpc" 1 (B lh + B 3h ) contributes ~ 10% of the total 
bispectrum, while at k ~ 10 h Mpc" 1 B%h contributes less than 1% and B\h less than 5%. This happens 
both for the bispectrum with Gaussian initial conditions and for the non-Gaussian component, AB^q. 
Thus, the following discussion will be valid in the squeezed limit, independently of the initial conditions. 

In particular, let us consider the squeezed (not necessarily isosceles) configuration fa <C &2 — ^3 = k. 
If fa is still in the linear regime, then we can safely set that the Fourier transform of the halo profile is 

p(m, z, fa) ~ m (50) 



in all bispectrum terms involving fa in Eqs. (|15h , (16) and (17). By making this substitution, the 



expressions for the 1-, 2- and 3-halo terms greatly simplify. Moreover, by using the conditions in Eq. (35) 



6 Note that the effect oc 1/k 2 due to non-Gaussian corrections in the linear bias does not blow up the amplitude 



of the matter power spectrum at large scales thanks to the constraint Eq. ( 49 1 , which applies at those scales 
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Figure 3: Non-Gaussian corrections to individual ingredients of the Halo Model for the matter power 
spectrum (left panel) and for the squeezed configurations of the matter bispectrum (right panel). In 
addition to the full Halo Model (continuous curve) we consider the results of allowing for a non-Gaussian 
correction only in each one of the ingredients of the Halo Model: the mass function (short- dashed curve), 
the halo density profile (dot-dashed), the linear halo bias (dotted) and, for the bispectrum in particular, 
the quadratic halo bias (medium- dashed) and the initial component to the tree- level bispectrum (long- 
dashed). The correction to the quadratic bias (upper right panel) is shown with the sign changed. All 
panels share the same label. 



and Eq. (37), together with Eq. (50), the Halo Model bispectrum contributions become, at leading order 
in ki/k -C 1, 



1 



B\h(k\,k, k) = -e 2 s (k, /nl), 
B 2h (h,k,k) = e[ bl] (k,h L )P L (k l ) 



B 3h (k h k,k) = 2 



13 
14 



+ 



4 ldlnP L 
7 ~ 2 dink 



with ^j^jr evaluated at k. The functions in these expressions are defined as 

1 



(51) 
(52) 

P L (h) P 2h (k) , (53) 



ef ] (k,h L ) 



P' 



dmn NG (m, z, / NL ) /3 l (m, z, fe, /nl) ^(m, z, /nl) , 



(54) 



where F(m, z, /nl) represents a generic function of mass and redshift. Thus, these functions are like an 
"average" of the function F, weighted by the mass function and the ith power of the Fourier transform 
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Figure 4: Left panel: functions ef 1 \k) (continuous curve), elf (A^/elj (0) (dashed curve), ef 2 \k) (dot- 
dashed curve) and e 2 m \k) (dotted curve) as a function of k, evaluated for Gaussian initial conditions at 
redshift z = 0. Right, top panel: non-Gaussian to Gaussian ratio for the functions e^\k) (continuous 
curve), e 2 bl \k) (dashed curve). Right, bottom panel: comparison between the function e 2 bl \k) / e 2 bl \o) for 
Gaussian initial conditions (dot-dashed curve) and non-Gaussian initial conditions (long dashed curve). 



of the density profile. The first two terms inside the bracket of Eq. (53) have been derived by taking the 
squeezed limit k\ <C k 2 ~ k 3 of Bq(Ux, k 2 , k 3 ) in Eq. i.e., 



B G (h, k 2 , k 3 ) ~ 2 [F 2 (ki, k 2 )P L (ki)P L (k2) + F 2 (k 1 ,k 3 )P L (k 1 )P L (k 3 )] 



13 
14 



+ 



4 ld\nP L 
7 ~ 2 dlnk 2 



(k!-k 2 ) 2 + 0(h/k 2 ) 



P L (h)P L (k 2 ), 



(55) 



where we have expanded k 3 = |ki + k2 1 for small k\/k 2 . Note that only B^, Eq. (53), depends on the 
angle between the long and the short mode. For isosceles configurations ki • k 2 ~ 
angular dependence drops. 



0(h/k 2 ) and the 



Eqs. (51 )-(53) also hold for Gaussian initial conditions. In this case the term 2 /nl/-^(^i> z ) inside the 
bracket drops and n^Q(m, z, /nl), p(m, z, k, /nl) and F(m, z,/nl) are replaced by their corresponding 
Gaussian expressions. Moreover, these equations can be also extended to other types of non-Gaussian 
initial conditions, provided that the corrections to n^c, p, b\ and b 2 are known. The derivation is 
analogous to the one given here for local non-Gaussianities. 

In Fig.[4](/e/i panel) we plot the functions involved in Eqs. (51 )-(53), evaluated for Gaussian initial 
conditions. As for small k, p(m, z, k) 

0. For k < 0.2/iMpc" 1 



i.e. e 



1 and e? 2j 



m, then 



and ef 2 ^ reduce to the constraints of Eqs. (j35h — (|3Th 



only deviates less than 0.1% from unity while 



is smaller than 0.001. From these observations we conclude that the approximations above, Eqs. (51) 



(53), are safely valid when k\ < 0.2/iMpc -1 . On the right panels of Fig.|4]we consider instead the effects 
of non-Gaussianity on the functions. In particular, the upper right panel shows the ratio between ef 1 ^ 
for /nl = 100 to the same quantities for Gaussian initial conditions. We notice that the effect on ef 1 ^ is 
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Figure 5: Comparison between the full Halo Model prediction of Eqs. (15)-(17) (thin, black curves) 
and the approximate expressions of Eqs. (51)-(53) (thick, blue curves) for squeezed configurations with 
k\ = 0.014 /iMpc -1 as a function of k\ = k% = k, at redshfit z = 0. Dotted, dashed and dot-dashed curves 
correspond to the 1-, 2- and 3-halo contributions while the continuous curve correspond to their sum. On 
the left panel we show the Gaussian prediction Bq while on the right panel we show the non-Gaussian 
component A1?ng = -Bng — Bq . The lower panels present the ratio between the approximations and the 
full expressions for each component. The shaded area marks value within a 10% discrepancy. The thin 
vertical line indicates the value of k\ , and therefore corresponds effectively to an equilateral bispectrum 
configuration. 



very small, below one percent over the range of scales considered, while for e 2 it is more significant, of 
the order of 5% at large k. The lower right panel compares ef 2 ^ for /nl = 100 and /nl = 0, illustrating 
that non-Gaussian initial conditions affect significantly this function at small scales. 

In Fig. [5] we show a comparison between the full Halo Model prediction (thin, black curves) and 
the approximate expressions presented above (thick, blue curves) for squeezed configurations with k\ = 
0.014 /iMpc -1 as a function of k2 = k% = k, at redshfit z = 0. The left panels present the matter bispec- 
trum with Gaussian initial conditions Bq, while the right panels present the non-Gaussian component, 
A£>ng = -Bng — Bq. The lower panels show the ratio between the approximate and the full expressions, 
distinguishing the 1-halo, 2-halo and 3-halo contributions. 

As expected, the approximate predictions describe well the exact evaluations for large values of 
k. More precisely, the difference between the two is below 10% for k > 0.02 /iMpc -1 with k\ = 
0.014 h Mpc -1 , that is, already for a ratio between the two sides given by k/k\ > 1.5. For kjk\ > 10 the 
error is at the percent level for Bq and even lower for A1?ng- However, limiting our attention to the 
approximation of the i?2/t term alone, we notice that one needs to consider at least k/k\ > 10 to reach 
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an accuracy below 10%. As expected, the approximation for B\h works better than 0.01% for all the 
momenta range of interest, both in the Gaussian as in the non-Gaussian case. 

To conclude, for k\ < 0.2/iMpc -1 and k > fei, the total bispectrum in the squeezed- limit, k\ <C 
k = k 2 — k3, with Gaussian (/nl = 0) or non-Gaussian (/nl 0) initial conditions, can be very well 
described by the expression 



1 



B(h, k,k) = -_ t [ p (k, / NL ) + \ 4 6lJ (k, /nl) + 2 P 2h (k, / NL 
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to the ~ 10% and ~ 1% levels for k/k\ > 1.5 and k/k\ > 10, respectively. Since at small scales 
(k > 0.2/iMpc _1 ) and for squeezed configurations B 2 h dominates over Bih and B$h, for k > 6/iMpc -1 
we find that the total bispectrum is well described (to better than 1% as compared to the full Halo 
Model) simply by 



1 



B(k\, k,k)~-_ e [ p (k, / NL ) + e[ bll (k, /nl) P L (h 



P 



(57) 



where we have included the term e 2 ™ (k, /nl)/p (coming from the approximation of Bih) to guarantee a 
accuracy better than 1%. Should we drop it and consider only the term e£ (k, /nl) -Pl(&i)> then the Bis- 
pectrum would be captured at a 5% level for k >QhMp C -\ Only for k > 30 /iMpc" 1 , ef l] (A:, / NL ) P L (ki) 
alone describes, to better than 1%, the full matter bispectrum in the squeezed configuration. 



4.3 Equilateral and generic triangles 

Figure [6] shows the same results as in Fig.[2]but for equilateral configurations, k\ = k 2 = k%. In particular, 
we plot B(k,k,k) as a function of k = k\ = k 2 = k% using the same line style as in Fig. [2j The first 
row shows the matter bispectrum for Gaussian initial conditions, together with the ratio between the 
measurements in N-body simulations and the model. The second and third rows show, respectively, 
the difference and the ratio between the non-Gaussian and Gaussian case. For z = the agreement is 
reasonably good, especially on small scales. However, for z = 1 the agreement worsens and for both 
redshift it does not remain within the 10% accuracy found in the squeezed limit. On the other hand, it 
is reasonable to expect that such large discrepancy is limited only to the mildly nonlinear regime, where 
the HM typically fails. 

The non-Gaussian components represent a correction of a few up to ~ 10%. An enhancement, 
shown as a bump in the ratio between the non-Gaussian and Gaussian bispectrum, is also displayed for 
equilateral configurations. Such local maximum is located at similar but slightly bigger scales w.r.t. the 
squeezed case. We find it at k w 0.65 h Mpc -1 for z = and at k ~ 1.0 h Mpc -1 for z = 1. The difference 
in the correction when comparing the bispectrum for z = versus the one for z = 1, is not as big as in 
the squeezed configuration, where there was a change from a 15% to 25% enhacement. Here we rather 
observe a change in the correction from 9% at z = to 11% at z = 1. 

The ratio of the difference -Bng — Bq between the N-body simulations and our predictions is inside 
the 10% agreement for z = 0, but largely outside for z = 1. We do not attempt here to explain why 
the Halo Model is less accurate for equilateral configurations when local non-Gaussianity is considered. 
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Figure 6: Equilateral configurations of the matter bispectrum, B(k, k, k), at redshift z = (left panels) 
and z = 1 (right panels). See text for explanation. 
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Figure 7: Relative effect of primordial non-Gaussianity shown in terms of the quantity 
^ng(^I) &2i kz)/BG(ki, ^2, k$) at z = 1 plotted as a function of the ratios k^/ki and k^/kx assuming 
a constant k\ = 3/iMpc — 1 . The function takes values on a triangle (represented by the shaded area on 
the bottom surface) where the lower corner corresponds to flattened triangles (k\ = 2k2 = 2ks), to top 
right corner to equilateral configurations (k\ = ki = k%) and the top left corner to squeezed configurations 
{ki <C k\ = ks). The left panel shows the ratio B^q/Bq computed in the three-level approximation in 
PT. The right panel shows instead the full HM calculation. 



However, since the ratio B^q/Bq agrees very well with the one measured in the numerically simulations, 
this inaccuracy is likely to be traced in the Halo Model itself, more than on the inclusion of non-Gaussian 
corrections. 

In order to provide a more complete description of the impact of primordial non-Gaussianity on all 
triangular configurations, in Fig. [7] we show the relative effect Sng(&i 5 ks)/BQ(k\, &2> ks) at 2 = 1 
as a function of the ratios ks/ki and A^/fci assuming a constant k\ = 3/iMpc . In order to avoid 
redundancy among equivalent configurations, the quantity plotted takes values on a triangle (represented 
by the shaded area on the bottom surface) where the lower corner corresponds to flattened triangles 
(k\ = 2k2 = 2^3), the top right corner to equilateral configurations (k\ = k<i = £3) and the top left 
corner to squeezed configurations [k^ *C k\ = £3). The left panel shows the ratio B^q/Bq computed in 
the tree-level approximation in PT. Even though perturbation theory is not applicable at these scales, 
this plot nevertheless provides an estimate of the contribution of the linear bispectrum Bo to the overall 
effect. Clearly, a large effect is present only in the squeezed limit. The right panel shows instead the 
full HM calculation. In this case we notice how, in addition to a 20% effect in the squeezed limit, all 
triangles of essentially any shape receive a correction of the order of about 7-8%. We also notice, as 
in the previous plots, a local maximum effect corresponding to wavenumbers of order 1/iMpc -1 and to 
flattened configurations. 
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5 Conclusions 



In this paper we have used the Halo Model approach to compute, over a wide range of scales including 
very small ones, the matter bispectrum for non-Gaussianity of the local type. We have considered non- 
Gaussian corrections which appear in the three fundamental ingredients of the Halo Model: the halo 
profile, the halo mass function and the bias functions. For the Gaussian case we have used the NFW 



(Navarro et al. 1997) halo profile and the Sheth & Tormen (ST,|Sheth fc Tormen 1999) expressions for 



the halo mass and bias functions. The non-Gaussian correction to the halo profile has been included by 
using Eq. (22), which is derived as a fit to the numerical simulations presented in Smith et al. (2011). 



For the correction to the halo mass function we have employed the expression proposed by Lo Verde 



et al. (2008) based on the Edgeworth expansion about the Press-Schechter mass function, with v = 5 c /a 
rescaled by a scaling parameter q = 0.91. Finally, for the bias functions we have assumed a simplified 



version of the expressions in Baldauf et al. (2011) and Sefusatti et al. (2011), as described in sec. 3.3 



Using these ingredients, we have computed the matter power spectrum and bispectrum, both in the 
Gaussian and non-Gaussian cases. For a local parameter /nl = 100, in Fedeli & Moscardini (2010) 
and Smith et al. (2011) it was noted that non-Gaussian corrections enhance the matter power spectrum 
on scales k ~ l/i _1 Mpc, inducing a bump in the ratio between the non-Gaussian and Gaussian power 
spectrum. As shown in Fig. [3j this bump is mainly due to the corrections to the halo mass function, 
and represents a 2 — 3% correction. For the same parameter /nl = 100, we have shown that the non- 
Gaussianity introduces a correction in the squeezed configuration k\ <C k<i ~ k% = k of the bispectrum, 
of about ~ 15% for redshift z = and ~ 25% for redshift z = 1, also at k ~ 1 /iMpc -1 . We thus confirm 
the presence of a similar bump also in the ratio between the non-Gaussian and Gaussian bispectra. We 
have also shown that the non- Gaussian corrections in the equilateral configuration k\ = ki = k% = k of 
the bispectrum, reach a maximum at k ~ 1 /iMpc -1 of about ~ 10%, both for redshift z = and z = 1. 
The corrections due to local non-Gaussianity in the equilateral configuration if the bispectrum are thus 
less sensitive to the redshift. 

We have also studied how the non-Gaussian corrections affect each of the halo model ingredients 
individually, which is shown in Fig. [3j This has allowed us to conclude that, contrary to the case of 
the matter power spectrum, it is the non-Gaussian correction in the linear bias which determines the 
shape and amplitude of the maximum distortion of the squeezed configuration of the Bispectrum, in the 
presence of local non-Gaussianity. Moreover, the corrections to the halo profiles are subdominant in the 
bispectrum over all the scales considered in this paper. Finally, we conclude that all corrections of the 
Halo Model ingredients might be important, since different correlators seems to be sensitive to different 
corrections of the ingredients. 

We have also compared our results to numerical simulations both for squeezed and equilateral configu- 
rations. In the squeezed limit the agreement is remarkably good, both in the Gaussian and non-Gaussian 
case. In particular, it is better than about 10% for z = 0, both for the total bispectrum and for its 
non-Gaussian correction. This strongly suggests that in the squeezed limit the Halo Model predictions 
for the matter bispectrum work reasonably well and they capture with very good accuracy non-Gaussian 
corrections of the local type. For equilateral configurations the agreement between the Halo Model and 
N-body simulations is still reasonably good for z = but worsens for z = 1. 

In real space, the squeezed limit represents the case where two points are close to each other, belonging 
to the same halo, and far from the third point, belonging to a different halo. In this limit we have checked 
that the bispectrum is correctly dominated by the 2-halo term. Moreover, on large scales the Fourier 
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transform of the halo profile in the long mode can be simply approximated as the mass of the halo. This 
allows a considerable simplification of the Halo Model expressions, which can be described by very simple 
formulas, Eq. (56) at mildly non-linear scales and Eq. (57) at highly non-linear scales. These represent 
master formulas for the matter bispectrum in the squeezed configuration independently of the nature of 
the initial conditions. 

Given the agreement that we have found between the Halo Model calculation and the simulations for 
the matter bispectrum, it seems promising to use the Halo Model to compute observable quantities such 
as the galaxy or lensing bispectra. We indeed provide a preliminar estimate of the galaxy bispectrum in 
the presence of local non-Gaussianity in Appendix [XJ 



Acknowledgements 

D. G. F. is supported by the Academy of Finland grant 131454 and also acknowledges support from a 
Marie Curie Early Stage Research Training Fellowship associated with the EU RTN 'UniverseNet', during 
his stay at CERN TH-Division when this project was initiated. A.R. is supported by the Swiss National 
Science Foundation (SNSF), project 'The non-Gaussian Universe" (project number: 200021140236). E.S. 
was supported in part by the EU Marie Curie Inter-European Fellowship. 



A The Galaxy Bispectrum 



The Halo Model framework can be extended to predict galaxy correlators by means of the Halo Occupation 
Distribution (HOD). In fact, the additional ingredient required is given by a prescription describing the 
way galaxies populate dark matter halos. In analogy with Eq. (|8|), we could write the galaxy power 



spectrum as the sum (Seljak, 2000, Scoccimarro et al. 2001) 



P g (k)=P lh>g (k) + P 2h , g (k), 
where now the two contributions are given by 

(N g (N g - l)) m p 2 (k,m) 
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P h (k,mi,m 2 ) , 



(58) 

(59) 
(60) 



with n g being the mean galaxy number density and where {N g ) m and (N g (N g — l)) m represent, respec- 
tively, the first and second moments of the probability P{N\m) for a halo of mass m to contain N galaxies. 
In simple terms, {N g ) m describes the mean number of galaxies in halos of mass m while {N g {N g — l)) m 
describes instead the mean number of galaxy couples in the same halo. Similarly, for the bispectrum we 
have 

B g (ki,k 2 , h) = B lh>g (k\, k 2 , k 3 ) + B 2h>g (ki, k 2 , k 3 ) + B 3htg (ki, k 2 , fc 3 ), (61) 
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(64) 



where, in addition, the 1-halo term now requires the third momenlj^jof the HOD, (N g (N g — l)(N g — 2)) m . 

A given HOD clearly depends on the specific galaxy population that one wants to describe. In the 
following, for illustration purposes we will consider the case of Luminous Red Galaxies (LRG), since this 
is today the best know galaxy type. It has been shown that a proper description of the LRGs HOD 
can be obtained by separately considering, for halos containing more than one galaxy, a central galaxy 
(i.e. a single galaxy placed at the halo center) and a population of satellite galaxies following a Poisson 



distribution (Kravtsov et al. 2004). In particular we will adopt the following parametrization for the 
mean 

(N g ) m = (N c ) m + (N s ) m , (65) 

with 
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where m m ; n represents the mass threshold above which we expect at least one, central galaxy while 
m\ is the typical mass above which a second galaxy is present. We assume the values of the three 
parameters m m ; n , m\ and a from Kulkarni et al. (2007), which obtain them from a fit to measurements 



of the LRGs three-point function in the Sloan Digital Sky Survey. Their best fit corresponds to m mm = 
7.66 x 10 13 h^ 1 M , m\ = 4.7 x 10 14 h~ l M and a = 1.4. In order to compute the power spectrum and 
bispectrum 1-halo terms we will need as well the second and third order moments. Since we assume the 
satellite HOD to correspond to a Poisson distribution, then we simply have 



(N s (N s -l)...(N s -j)), 



so that 
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(68) 
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7 Following the review of Cooray & Sheth ( 2002| ) one should notice that: 1) in principle the spatial distribution 
of galaxies inside the halo does not necessarily follow the dark matter distribution, but it is indeed a good approx- 
imation; 2) if an halo contains a single galaxy then this is expected to sit at the center of the halo itself, and the 
double power of the profile should be replace by a single one. These details are not essential for our purposes. 
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In these calculations, the problem represented by the constant value assumed by the 2-halo and 1-halo 
contributions at small k is more severe than in the matter case. In fact, for both squeezed or equilateral 
configurations such contributions are not negligible at large scales. A possible solution to this problem 
has been proposed by assuming "compenpensated profiles" for the halo matter ditribution (see, e.g., 



section 4.4 in Cooray Sz Sheth, 2002). In our illustrative calculations we adopt the simple prescription 
consisting in replacing the NFW profile p(m, k) with the difference p(m, k) — mWn(k), Wn(k) being the 
top-hat window function of radius R = [3m/ {A-np)] 1 ^ . We perform such replacement only for the profiles 
with /c2 ~ k% = k in the expression for the 1-halo contribution, Eq. (62), and for those in the first integral 



in the 2-halo contribution in Eq. (63). While this is clearly not a well justified procedure, our intent is to 



avoid a largely incorrect behaviour at large scales and focus on the small scale non-Gaussian correction, 
which is not affected by this choices. 

In Fig.[8]we show the results for squeezed configurations of the galaxy bispectrum B g (k\, k, k), defined 
by ki = 0.014 /iMpc" 1 fixed, plotted function of k 2 = k% = k. Similarly to the matter case, left 
panels show the results at redshift z = 0, while right panels at z = 1. In the first row we present the 
galaxy bispectrum for Gaussian initial conditions. The middle panels show the component due to non- 
Gaussian initial conditions with a local /nl = 100, that is the difference between the predictions with 
and without non-Gaussianity, AB^ G = -B^ g (/nl = 100) — -B g (/nl = 0) while the lower panels show 
the non-Gaussian to Gaussian ratio. In all plots, the dotted, dashed and dot-dashed curves correspond, 
respectively, to the 1-halo, 2-halo and 3-halo contributions, while the continuous line corresponds to their 
sum. 

Again the dominant term in this squeezed limit is given by B2h(ki,k,k). It is interesting to notice 
how, despite presenting as well a maximum around k — 1 /iMpc -1 at z = 0, unlike in the matter case, the 
ratio between the non-Gaussian galaxy bispectrum to the Gaussian one, shows a more constant behaviour 
as k grows. The non-Gaussian effect is, in fact, of the order of 20% (30%) over a large range of scales 
at redshift z = (z = 1). This is a significant correction that might be taken into account in studies of 
galaxies bias at small scales. 

As in the matter case, for the galaxy bispectrum as well we can find simplified expressions for squeezed 



triangular configurations. The constraints of Eq. (|35[) and (37), will be now modified into 
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where b\ tg and 62, g represent the linear and quadratic bias of the galaxy population. Taking into account 
these relations we can obtain now a squeezed approximation for the galaxy bispectrum, using, as we did 
in the matter case, the fact that p{k) —> m when k — > 0. For k\ <C k we find 
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Figure 8: Squeezed configurations of the galaxy bispectrum, B g (k±,k,k), with k\ <C k<i = = k, at 
redshift z = (left panels) and z = 1 (right panels). See text for explanation. 
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Figure 9: Comparison between the full Halo Model prediction of the galaxy bispectrum of Eqs. (|62|), (|63|) 



and (64) (thin, black curves) and the approximate expressions of Eqs. (73), (74) and (73) (thick, blue 
curves) for squeezed configurations with ki = 0.014 /iMpc^ 1 as a function of kx = kg = k, at redshfit 
2 = 0. Dotted, dashed and dot-dashed curves correspond to the 1-, 2- and 3-halo contributions while 
the continuous curve corresponds to their sum. On the left panel we show the Gaussian prediction -B s ,g 
while on the right panel we show the non-Gaussian component ABg^c = B g ^c — B g c- The lower 
panels present the ratio between the approximations and the full expressions for each component. The 
shaded area marks value within a 10% discrepancy. The thin vertical line indicates the value of k%, and 
therefore corresponds effectively to an equilateral bispectrum configuration. 



where the e„ ■ functions read now 



^- / dmn(m, / NL ) 
nl J 



p(k,m,f NL ) 



m 



F(m, /nl) 



(76) 



In Fig. [9]we compare the simplified expressions Eqs. ( 73 )— ( 75 ) versus the exact Halo Model evaluation, 

+100. The left panels consider the case of Gaussian 



Eqs. (|62|)-(|64|), in the case of local NG with 



initial conditions while the right panels show the non-Gaussian difference, as done in Fig. [5] for the 
matter bispectrum. It is evident that in order to achieve an agreement better than 10% it is necessary, 
in all cases and for every individual term, to consider k > 0.5/iMpc -1 for ki = 0.014 h Mpc -1 , that is 
a ratio between the two sides of k/k\ > 40. For k > 1 foMpc -1 (k/ki > 100), the simplified expressions 
reproduce the full ones at better than 1% level. 
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B Details on the mass integrals 



In this work all mass integrals in the HM expressions are performed numerically. An appropriate upper 
bound to such integration is easy to find, since the halo mass function decays exponentially fast. On 
the other hand it is clearly not feasible to set to zero the integral's lower bound. We are then forced 
to cut the integration at some minimum value M m [ Q and compensate for the distribution low-mass tail 
removed for M < M m [ n . In this Appendix we explain our prescription for this procedure and show that, 
for sufficiently small values of M m ; n , our results are insensitive to such cut-off scale, as it should be. This 
is particularly important for the constraint equations, Eqs. (35)-(37). To keep general the discussion 
we will then discuss the non-Gaussian case, while the Gaussian case can simply be obtained by setting 
/nl = 0. 

Let us begin by the conservation of mass contraint equation, 

1 



<dMMn(M,z;f NL ) = l. (77) 

P J 

Since we will only integrate over a finite range [M m ; n , M max ], we can compensate for the removed low 
mass tail distribution at M < M m j n , by introducing a constant Co, such that 

If 1 f M ™* M ■ 

- dM Mn(M,z;f NL ) = - dM M n(M, z; / NL ) + C ^^ = 1. (78) 

P J P JM min P 

This is equivalent to consider the mass function, for M ^ M m ; n , as a Dirac-delta like n(M,z; /nl) = 
Cq(z, /nl; M m j n )5D(M — M m ; n ). The constant Co can be easily obtained from the above equation, as 

C (z, / NL ; M min ) = [1 - D (z, / NL ; M min )]-^- , (79) 

-'"min 

D (zJ NL -M miQ ) = - dM M n(M, z; /nl) • (80) 

P J Mmin 

For any choice of M m j n there will be always a constant Co, given in terms of Do, that will just act 
in such a manner that the mass contraint will be attained like in Eq. (78). The constant Dq measures 
the deviation from unity when the constraint equation is evaluated by integrating within [M m i n , M max ]. 
Theorefore, for any integration over the mass in the Halo Model expressions involving the product of 
the mass function n(M) times some other function F{M) (except the bias functions, see below), we will 
simply apply the prescription 

ldMn{M,z;f NL )F(M;...) = I dM M n(M, z; /nl) F(M; ...) 

JM min 

+ C (z,h L ;M min )-^F(M min ;...). (81) 

-'"min 

We have found that the predictions from the Halo Model are insensitive to decreasing further the lower 
bounds below M m - m = 10 6 /i -1 M , so we have taken this value as the minimum mass cut-off. 
Other constraint equations, involving the bias functions, are given by 
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where A6^ G are the scale-independent (SI) corrections to bi and Abf^ G the scale-dependent (SD) ones. 
Again, we enforce these constraints by adding extra terms as 

^ jdmmn^c (b ijg + A6^ G ) 

= \ [ dmmn^G (b ijg + Ab^ G ) + cf l \m, z, / NL ; M mm )^^ = 5a , (84) 



3/ dmmn^G Abfl G + c} SD) (m, z, /nl, fc;M min )^^ = , Vz ^1, (85) 

P JM min P 



where the new constants read 



Cf l) (z, / NL ; M min ) = [5 a - Df l) (z, /nl j M min )]^ , (86) 
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Integrating over the mass the product of the mass function n(M) with any bias function (or its non- 
Gaussian corrections) times some other function F(M), we apply an analogous prescription like 



dMn NG (\ g + A6g G ) F(M; . 

J 

J A 



na dMMn NG (b hg + Ab^ G ) F(M,...) + C'f I) (^/ NL ;M min ) T /- F( ^ mui :..,). (90) 



M min 

(b) 



dMn NG A6^ G F(M;...) 



/ 



dMMn NG Ab\% F(M, ...) + cf D) (z, /nl, fc; M min )-J— F(M min ; ...). (91) 



We have checked that the predictions from the Halo Model where the addition of these extra constants 
are involved, are again insensitive to decreasing further the lower bound below M m \ n = 10 6 h~ x M . Thus 
we have choosen that value as the minimum mass cut-off for any mass integral all the Halo Model terms. 
Finally note that after introducing Co, the additional assumption of Cj(SI) and Cj(SD) was necessary 
and not redundant, since the latter are not directly related to Cq. Should any other constraint equation 
in the form of a mass integral over the mass function or the bias functions times other (mass dependent) 
functions exist, then another 'compensating' constant should be defined following the same philosophy. 
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